Software Implementation of Cellular Automata 
Based Random Number Generator 
J. Barry SHACKLEFORD et al. 
HP Docket No. 10019023-1 
1/14 



LUT ADDRESSES 
8 4 2 1 



0 
0 
0 
0 
0 
0 
0 
0 



0 
0 
0 
0 
1 
1 
1 
1 
0 
0 
0 
0 
1 
1 
1 
1 



0 
0 
1 
1 
0 
0 
1 
1 
0 
0 
1 
1 
0 
0 
1 
1 



0 
1 
0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 



r 

0 

1 
1 
1 

0 
0 

1 

0 

1 

1 

0 

1 
1 

0 
0 
0 



1 



0001101101001110 



==> 1 B4E 16 
Z=> 6 9 9 0 10 
Z=> CA06990 



FIG 1A 




=> 




FIG 1B 



Software Implementation of Cellular Automata 
Based Random Number Generator 
J. Barry SHACKLEFORD et al. 
HP Docket No. 10019023-1 
2/14 



i+d f 



i+d 4 




i+d 2 := ^^ ( d s> d 4> d 2> d i) 
i+d 1 

i 



FIG 1C 



(-1,0,1,2) 




O O O 



FIG 1D 



Software Implementation of Cellular Automata 
Based Random Number Generator 
J. Barry SHACKLEFORD et al. 
HP Docket No. 10019023-1 
3/14 



2QQ 



DETERMINE 
EMULATION 
PARAMETERS 

210 




f 


INITIALIZE 




220. 






RUN SIMULATION IN 
PARALLEL 




230 




r 


OUTPUT RESULT 




240 



FIG 2 



Software implementation Of Cellular Automata 
Based Random Number Generator 
J. Barry SHACKLEFORD et al. 
HP Docket No. 10019023-1 
4/14 



#include <stdio.h> 
#include <sys/types • h> 
#include <sys/uio.h> 
#include <unistd.h> 
#include <fcntl.h> 
#include <string.h> 



300 



#ifndef BUFSIZ 
#define BUFSIZ 4096 
#endif /* BUFSIZ */ 
#define WSIZE 32 

#define scmp(a,b) (strcmp ( (a) , (b) ) == 0) 
#define rotateLef t (a, b) ( (a<<b) | (a>> (WSIZE-b) ) ) 



/* Logic for TruthTable 51510 */ 

#define f <a8 , a4 , a2 , al) ( (a8 & a4 & a2 & al ) | (a8 & a4 & a2 & -al ) | (a8 & ~a4 
& a2 & al ) | (a8 & -a4 & -a2 & -al ) | (~a8 & a4 & -a2 & al ) | (~a8 & a4 & ~a2 & 
-al ) | (-a8 & -a4 & a2 & -al ) | (~a8 & ~a4 & -a2 & al ) ) 



static unsigned long stateEven, stateOdd ; 
unsigned long * returns tate = &stateEven ; 

char *outFile = "bsrand.32" ; 
unsigned long initUpper = 0x80000000 ; 
unsigned long initLower = 0x0 00 0 0000 ; 
long initRun =64 ; 
long numData = 3000000 ; 

void bsinit (unsigned long upper, unsigned long lower) 

{ 

int i ; 

stateEven = 0 ; 
stateOdd = 0 ; 

for ( i = 0 ; i < 32 ; i += 2 ) { 
stateEven <<= 1 ; 
stateOdd <<= 1 ; 

stateEven |= (upper >> (31 - i) ) & 0x01 ; 
stateOdd |= (upper >> (30 - i) ) & 0x01 ; 

} 

for ( i = 0 ; i < 32 ; i += 2 ) { 
stateEven <<= 1 ; 
stateOdd <<= 1 ; 

stateEven |= (lower >> (31 - i) ) & 0x01 ; 
stateOdd |= (lower >> (30 - i) ) & 0x01 ; 

} 

} 

/* { -7, 0, 11, 17 } */ 
unsigned long bsrandO 
{ 

unsigned long newStateEven ; 
newStateEven = f ( 

rotateLeft (stateOdd , 28), 

stateEven, 

rotateLeft (stateOdd , 5), 
rotateLeft (stateOdd , 8) 

) ; 

stateOdd = f ( 

rotateLeft (stateEven, 29), 
stateOdd, 

rotateLeft (stateEven, 6) , 
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rotateLef t (stateEven, 9) 

) ; 

st ate Even = newStateEven ; 

return { *returnState ) ; 300 

} 

void usage (char *name) 

{ 

fprintf (stderr, "usage: %s name) ; 

f pr int f (stderr, " [-o outFileName] [-seed upper lower] ") ; 

fprintf (stderr, "[-even | -odd] [-help] " ) ; 

fprintf (stderr, " [-num numData] [-skip numlnitialRun] \n" ) ; 

} 

void getOpt(int argc # char *argv[]) 

{ 

while (argc > 1) { 

if ((argc > 2 ) &&scmp (argv [1] , n -o")) { 

outFile = argv [2] ; 

argc -= 2 ; argv +=2 ; 
} else if ((argc > 2) &&scmp (argv [1] , "-num")) { 

sscanf (argv [2] , "%ul" , &numData) ; 

argc --2 ; argv += 2 ; 
} else if ((argc > 2) &&scmp (argv [1] , "-skip")) { 

sscanf (argv [2] , "%ul" , &initRun) ; 

argc -= 2 ; argv += 2 ; 
} else if ({argc > 3) &&scmp (argv [1] , "-seed")) { 

sscanf (argv [2] , "%ul", SdnitUpper) ; 

sscanf (argv [3] , "%ul" , fcinitLower) ; 

argc - = 3 ; argv += 3 ; 
} else if (scmp (argv [1] , "-odd")) { 

returnState = &stateOdd ; 

argc-- ; argv++ ; 
} else if (scmp (argv [1] , "-even")) { 

returnState = fcstateEven ; 

argc-- ; argv++ ; 
} else { 

usage ( argv[0] ) ; 

exit ( 1 ) ; 

} 

} 

} 

int main(int argc, char *argv[]) 

{ 

int fd ; 

unsigned long num, i, j, unit ; 

unsigned long *buf ; 

unit = BUFSIZ / sizeof ( num ) ; 

buf = (unsigned long* ) malloc ( BUFSIZ ) ; 

getOpt ( argc , argv ) ; 

fd = open( outFile, 0_WR0NLY | 0_CREAT, 0666 ) ; 
bsinit ( initUpper, initLower ) ; 
for ( i = 0 ; i < initRun ; i++ ) { 
bsrandO ; 

} 

for ( i = 0, j = 0 ;i< numData ; i++ ) { 
buf[j++] = bsrandO ; 
if ( j == unit ) { 

write (fd, buf, BUFSIZ) ; 
j * 0 ; 
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} 

} 300 

close ( fd ) ; 
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/******************************************** r»r\r\ 

,^ ^ 600 

* mkbsrand.c 

* 

* make Barry Shacklef ord' s rand generator program. 

* $Header: /home /mot o /work /FPGA2 002 /DIEHARD/ src/MKBSRAND/RCS/mkbsrand . c , v 
1.2 2001/09/07 21:26:47 moto Exp moto $ 

* 

* Created: 2-Jul-2001 

* Modified: $Date: 2001/09/07 21:26:47 $ 

* Author: Motoo Tanaka 

* Copyright (c) Hewlett Packard Company 
★*******************************************/ 

#include <stdio.h> 

unsigned int truthTable = 0 ; 

int a [4] , d[4] ; 

int isFirst = 1 ; 

FILE *ofp ; 

FILE *efopen(char *filename, char *mode) 
{ 

FILE *fp ; 

fp = fopen( filename, mode ) ; 
if ( fp == NULL ) { 

fprintf (stderr, "Can not open %s\n" , filename) ; 

exit ( 2 ) ; 

} 

return ( f p ) ; 

} 

void usage ( char *name ) 

fprintf (stderr, "usage: %s truthTable dl d2 d3 d4 { outfile }\n", 
name ) ; 

} 

int script ( char *string ) 

{ 

return ( fprintf (of p, string) ) ; 

} 

void GetOpt(int argc, char *argv[]) 
{ 

int i ; 

ofp = stdout ; 
if (argc <= 5) { 

usage ( argv[0] ) ; 

exit ( 0 ) ; 

} 

if (argc > 5) { 

sscanf (argv[l] , "%lu", ^truthTable) ; 
sscanf (argv[2] , "%d», &a[0]) ; 
sscanf (argv [3] , "%d", &a[l]) ; 
sscanf (argv [4] , "%d" , &a[2]) ; 
sscanf (argv [5] , "%d", &a[3]) ; 
for (i = 0 ; i < 4 ; i++ ) { 

d[i] = a[i] ; 

if ( d[i] < 0 ) { 

d[i] += ( (-d[i] /64) + 1) * 64 ; 

} 

d[ i ] %= 64 ; 
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} 

If (argc > 6) { 600 
ofp = ef open (argv [6] , "w") ; 



void printHeadO 

{ 

script ( "#include <stdio.h>\n 1 ' ) ; 
script ( "#include <sys/ types .h>\n" ) ; 
script ( "#include <sys/uio .h>\n" } ; 
script ("# include <unistd . h>\n" ) ; 
script ( "#include <f cntl .h>\n" ) ; 
script ( "#include <string .h>\n n ) ; 
script ("\n") ; 

script ("#ifndef BUFSIZ\n") ; 
script ("#de£ine BUFSIZ 4096\n") ; 
script ("#endif /* BUFSIZ */\n") ; 
script ("#define WSIZE 32\n") ; 

script ( "#define scmp(a f b) ( strcmp ( (a) , (b) ) — 0)\n") ; 
script ("ttdefine rotateLef t (a, b) ( (a«b) | (a» (WSIZE-b) ) ) \n» ) ; 
script ("\n") ; 

} 

void printMiddle ( ) 

{ 

char string! 2 56 ] ; 
script ( n \n B ) ; 

script ( "static unsigned long stateEven, stateOdd ; \n" ) ; 
script ( "unsigned long *returnState = &stateEven ;\n") ; 
script ("\n") ? 

script ("char *outFile = V'bsrand. 32\" ;\n") ; 
script ( "unsigned long initUpper = 0x80000000 ;\n") ; 
script ("unsigned long initLower = 0x00000000 ;\n") ; 
script ("long initRun = 64 ;\n") ; 
script ("long numData = 3000000 ;\n\n" ) ; 
script ("void bsinit (unsigned long upper, "} ; 
script ( "unsigned long lower) \n{ \n" ) ; 
script ("\ tint i ;\n") ; 

script ("\tstateEven = 0 ;\n") ; 

script ("\tstateOdd = 0 ;\n") ; 
script ("\tf or ( i = 0 ; i < 32 ; i += 2 ) {\n") ; 
script ("\t\tstateEven <<= 1 ;\n") ; 
script ("\t\tstateOdd «= 1 ;\n") ; 

script ("\t\tstateEven |= (upper >> (31 - i) ) & 0x01 ;\n") ; 
script ("\t\tstateOdd |= (upper » (3 0 - i)) & 0x01 ;\n") ; 
script ("\t}\n") ; 

script ("\tf or ( i = 0 ; i < 32 ; i += 2 ) {\n") ; 
script ( "\t\tstateEven «= 1 ;\n") ; 
script ("\t\tstateOdd <<= 1 ;\n") ; 

script ("\t\tstateEven |= (lower >> (31 - i) ) & 0x01 ;\n") ; 
script ("\t\tstateOdd |= (lower >> (30 - i) ) & 0x01 ;\n") ; 
script ("\t}\n") ; 
script ("}\n\n") ; 

sprintf (string, "/* { %d, %d, %d, %d } */\n", a[0], a[l], a[2j, a[3]) 
script ( string ) ; 

script ( "unsigned long bsrand () \n" ) ; 
script C'{\n») ; 

script ("\tunsigned long newStateEven ;\n") ; 
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void printTailO 

{ 

script ( 11 \tstateEven = newStateEven ;\n" ) 
script ("\treturn{ *returnState ) ; \n") 
script ("}\n\n") ; 



HIS- 



600 



script ( "void usage (char *name) \n{\n" ) ; 

script ("\tfprintf (stderr, V'usage: %%s \", name) ;\n") ; 
script ("\tf print f (stderr, \" [-o outFileName] ") ; 
script (" [-seed upper lower] \") ;\n") ; 

script ("\tfprintf (stderr, \" [-even | -odd] [-help] \") ;\n") 
script ( "\tfprintf (stderr, \" [-num numData] " ) ; 
script (" [- skip numlnitialRun} \\n\ " ) ;\n") ; 
script ("}\n\n") ; 



script ( "void getOpt(int argc, char 
script ("\twhile (argc > 1) {\n") ; 



*argv[])\n{\n") 



script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 
script 



("\t\tif ((argc > 2) &&scmp (argv[l] , \»-o\ ,f )) {\n») ; 
("\t\t\toutFile = argv[2] ;\n") ; 
("\t\t\targc -= 2 ; argv += 2 ;\n") ; 

("\t\t} else if ((argc > 2 ) &&scmp (argv [1] , \"-num\")) {\n") , 
("\t\t\tsscanf (argv[2] , \"%%ul\", &numData) ;\n M ) ; 
( r, \t\t\targc -= 2 ; argv += 2 ;\n") ; 

("\t\t} else if ((argc > 2) &&scmp (argv [1] , V'-skipV')) {\n") 
("\t\t\tsscanf (argv[2] , \«%%ul\ n , fcinitRun) ;\n») ; 
("\t\t\targc -= 2 ; argv += 2 ;\n ,r ) ; 

(»\t\t} else if ((argc > 3) &&scmp (argv [1] , \"-seed\") } {\n") 

("\t\t\tsscanf (argv[2] , \ n %%ul\", ScinitUpper) ;\n") ; 

("\t\t\tsscanf (argv[3] , \ n %%ul\ f, f fcinitLower) ;\n") ; 

("\t\t\targc -= 3 ; argv += 3 ;\n") ; 

("\t\t} else if (setup (argv [1] , \»-odd\")) {\n n ) ; 

("\t\t\treturnState = &stateOdd ;\n") ; 

("\t\t\targc-- ; argv++ ;\n") ; 

("\t\t} else if (sctnp(argv [1] , \"-even\")) {\n") ; 

("\t\t\treturnState = &stateEven ;\n") ; 

("\t\t\targc-- ; argv++ ;\n") ; 

(»\t\t} else {\n») ; 

("\t\t\tusage( argv[0] ) ;\n") ; 

("\t\t\texit { 1 ) ;\n") ; 

(»\t\t}\n\t}\n}\n\n») ; 



script ("int main(int argc, char *argv [] ) \n" ) ; 
script ("{\n") ; 

script ("\tint fd ;\n") ; 
script ( "\tunsigned long num, i, j , unit ;\n") ,- 

script (" \tunsigned long *buf ;\n") ; 

script ("\tunit = BUFSIZ / sizeof( num ) ;\n") ; 

script ("\tbuf = (unsigned long* ) malloc ( BUFSIZ 
script ("\tgetOpt ( argc, argv ) ;\n") ; 

script ("\tfd = open( outFile, OJtfRONLY | O^CREAT , 

script ( "\tbsinit ( initUpper, initLower ) ;\n") ; 

script ("\tf or { i = 0 ; i < initRun ; i++ ) (\n") 

script ("\t\tbsrand() ;\n") ; 

script ("\t}\n") ; 

script ("\tf or ( i = 0, j = 0 ; i < numData ; i++ 
script ("\t\tbuf = bsrandO ;\n") ; 

script ("\t\t if ( j == unit ) {\n") ; 
script ("\t\t\twrite(fd, buf, BUFSIZ) ;\n") ; 
script ("\t\t\tj = 0 ;\n") ; 



) ;\n«) 



0666 ) ;\n") ; 
7 

) {\n«) ; 
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script ("\t\t}\n") ; 

script ("\t}\n") ; 

script ("\tclose ( fd ) ;\n") ; 
script ("}\n") ; 6 qq 



void cleanup () 

{ 

f close ( of p ) ; 

} 

* displace. c 

* Created: 2-Jul-2001 

* Author: Motoo Tanaka 

* Copyright (c) Hewlett Packard Company 



void nextEven ( int d ) 

□ { if ( d ) { 

,n if ( (d / 2)%32 ) { 

;S fprintf (ofp, "rotateLeft (%s, %2d) " , 

^1; (d % 2) ? "stateOdd " : "stateEven", d / 2 ) ; 

M } else { 

t fprintf (ofp, "%s», 

^ (d % 2) ? "stateOdd " : "stateEven" ) ; 

NJ } 

SO. } else { 

y fprintf (ofp, "stateEven") ; 

} 

□ } 

Ki void next Odd { int d ) 

° { if ( « . ( 

^ if ( ((d + 1) / 2) % 32 ) { 

fprintf (ofp, "rotateLef t (%s, %2d) " , 

(d % 2) ? "stateEven" : "stateOdd ", (d + 1) / 2 ) ; 
} else { 

fprintf (ofp, "%s" , 

(d % 2) ? "stateEven" : "stateOdd ", (d + 1) / 2 ) ; 

} 

} else { 

fprintf (ofp, "stateOdd") ; 

} 

} 

void nextState( int d[] ) 

{ 

int i ; 

fprintf (ofp, " \tnewStateEven = f (\n") ; 
for ( i = 0 ; i < 4 ; i++ ) { 

fprintf (ofp, "\t\t") ; 

nextEven ( d[i] ) ; 

if ( i 1= 3 ) { 

fprintf (ofp, ",") ; llKJI 

) 

fprintf (ofp, "\n") ; 

} 
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fprintf (ofp, "\t) ;\n») ; 

fprintf (ofp, "\tstateOdd = f (\n n ) ; 
for { i = 0 ; i < 4 ; i++ ) { 

fprintf (ofp, "\t\t") ; 

next Odd ( d[i] ) ; 

if ( i \= 3 ) { 

fprintf(ofp, ",") ; 

} 

fprintf (ofp, «\n") ; 

} 

fprintf (ofp, "\t) ;\n") ; 



* TruthTable to Logic 
* 

* Created: 2-Jul-2001 

* Author: Motoo Tanaka 

* Copyright (c) Hewlett Packard Company 
********************************************/ 

void i2b{ unsigned int data ) 

{ 

int mask ; 

for ( mask = 0x8000 ; mask > 0 ; mask >>= 1 ) { 
if (data & mask) { 

fprintf (ofp, "1") ; 
} else { 

fprintf (ofp, "0") ; 

} 

} 

} 



void emit ( int data ) 

{ 

int mask ; 

if ( isFirst ) { 

isFirst = 0 ; 
} else { 

fprintf ( of p x "|" ) ; 

} 

fprintf (ofp, « (») ; 

for ( mask = 0x8 ; mask > 0 ; mask >>= 1 ) { 
if ( mask < 0x8 ) { 

fprintf (ofp, "& ») ; 

} 

if ((data & mask)==0) { 

fprintf (ofp, "~") ; 

} 

fprintf (ofp, "a%X " , mask) ; 

} 

fprintf (ofp, ")") ; 

} 

void t21( unsigned int truthTable ) 
unsigned int mask ; 

FIG 

int i ; 
isFirst = 1 ; 
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fprintf (ofp, "/* Logic for TruthTable %u */\n" / truthTable) ; 
fprintf(ofp, "#define f (a8, a4, a2, al) (") ; 

for ( i = 15 , mask = 0x8000 ; mask > 0 ; mask >>= 1, i- - ) { 
if ( truthTable &. mask ) { 
emit ( i ) ; 

} } 600 

fprintf {ofp, ")\n") ; 

} 

void dolt ( ) 

{ 

printHeadO ; 
t21( truthTable ) ; 
printMiddle ( ) ; 
nextState { d ) ; 
printTail ( ) ; 

} 

main(int argc, char *argv[]) 

{ 

GetOpMargc, argv) ; 
doIt() ; 
cleanup () ; 
return { 0 ) ; 

} 
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